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The critical behaviour of helimagnets and stacked triangular antiferromagnets 
is analyzed in (4 — e) dimensions within three-loop approximation. Numerical es- 
timates for marginal values of the order parameter dimensionality N obtained by 
resummation of corresponding e-expansions rule out the possibility of continuous 
chiral transitions in magnets with Heisenberg or planar spins. 
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Chiral phase transitions in helical magnets and stacked triangular antiferromagnets with 
Heisenberg or XF-like spins as well as in some other systems attract much attention during 
the last decade |[]-||. Special interest to these transitions demonstrated both by theorists 
and experimentalists has been given rise by the conjecture that they belong to new 
universality class which is characterized by critical exponents differing markedly from those 
of 3D 0(n)-symmetric model with relevant values of n (2, 4, and 6). This conjecture orig- 
inates from the renormalizat ion-group (RG) analysis of corresponding (4 — e) -dimensional 
model performed within the lowest (one- and two-loop, according to quantum field theory 
language) orders in e. Some results given by the ^--expansion and Monte Carlo simulations 
were also considered as favoring the abovementioned idea 

In this Letter, an attempt is made to clear up whether the conclusion about an existence 
of new universality class for 3D chiral systems with Heisenberg or planar spins survives when 
higher-order terms in e-expansion are taken into account. Below we calculate e-expansions 
for quantities of interest up to the three-loop order. To obtain numerical estimates relevant 
to real 3D magnets we apply resummation procedures to these series before setting e = 1. 
Such a machinery proved to give good results for plenty of phase transition models. It is 
believed to be powerful enough to yield reasonable predictions in our case as well. 

The Landau-Wilson Hamiltonian describing systems under consideration may be written 
down in the form (see, e.g. Ref. ||): 

2 J x [ m ^ a ^ + V( ^ v< ^ + y¥Wa</W/3 + y<^a<A*</wJ > I 1 ) 

where <p a is a complex vector order parameter field, a, (3 — 1, 2, . . . , N, a bare mass squared 
ml being proportional to the deviation from the mean-field transition point. This model 
undergoes chiral phase transitions if wo > 0. In the opposite case, the transitions into 
somewhat trivial (linearly polarized or unfrustrated) ordered states take place. 

In the critical region, where fluctuations are strong and the system behavior is governed 



by the RG equations, the model Eq. ([[]) can demonstrate four different regimes of RG flow 
depending on N . Correspondingly, three critical (marginal) values of N exist separating 
these regimes from each other. If N < N& the RG equations possess three nontrivial 
fixed points (FP's) with the 0(2iV)-symmetric point being stable. When N exceeds iV cl 
the Heisenberg FP loses its stability but the other, anisotropic FP with coordinate w < 
acquires it. This point "annihilates" with another, saddle anisotropic FP when N approaches 
N C 2, and there is only one nontrivial FP in the domain N C 2 < N < N c ^. It is 0(2N)- 
symmetric and unstable. At last, when N increases further and crosses over the value iV C 3 
the creation of two new anisotropic FP's with w > takes place. One of them is stable and 
describes the chiral critical behavior. Hence, to answer the question about the relevance of 
the chiral FP to the critical thermodynamics of real helical magnets and stacked triangular 
antiferromagnets, one has to estimate and compare the number obtained with physical 
values N = 2 and iV = 3. 

Marginal values of N may be found analyzing RG /3-functions. We calculate these 
functions for the model Eq. ([!]) within three-loop approximation in (4 — e) dimensions using 
the minimal subtraction scheme (corresponding expansions are too lengthy and not presented 
here). The e-expansion for resulting from the /3-functions obtained is as follows: 

r ( 14v/6\ T137 91^6 /13 47 v / 6\ >/oN l , 
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= 21.80 - 23.43e + 7.088e 2 , (2) 

where ((x) is the Riemann ^-function, ((3) = 1.20206. The constant and linear terms in 
Eq. (H) coincide with those presented by H. Kawamura while the second-order one is 
essentially new. 

Such expansions are known to be asymptotic and physical information may be extracted 
from them provided some resummation method is applied. The Borel transformation com- 
bined with proper procedure of analytical continuation of the Borel transform usually plays 
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a role of this method leading to precise numerical estimates in cases of long enough orig- 



inal series [10,11 1. To perform the analytical continuation the Pade approximant of [L/l] 
type may be used which is known to provide rather good results for various Landau- Wilson 
models (see, e.g. Refs. p fi2yT3|[ ). The Pade-Borel summation of the expansion (fj) gives: 
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where a, b, and c are coefficients before e°, e 1 , and e 2 in Eq. (fj), respectively, Ei(x) being 
the exponential integral. Setting e = 1, we obtain from Eq. (|3]) 

iV c3 = 3.39 . (4) 

Making use of the Pade approximant [1/1] itself gives N c % = 3.81 while direct summation of 
the expansion Eq. (fj), being rather crude procedure, results in N c3 = 5.46. 

All these numbers, although considerably scattered, are nevertheless greater than 3. 
Hence, helical magnets and stacked triangular antiferromagnets with Heisenberg and XY- 
like spins are not seen to demonstrate new, chiral critical behavior. Instead, they should 
approach helically ordered state or frustrated antiferromagnetic phase only via first-order 
phase transitions. 

On the other hand, the difference between the number Eq. (Q) and iV = 3 is not so 
large. Moreover, numerical estimates for N c3 were obtained from the theory having no small 
parameter in the limit e = 1. How close to the precise value of iV C 3 they may be? 

To clear up this point let us compare Eq. (|4]) and its Pade counterpart with analogous 
estimate given by the RG analysis in three dimensions. Three-loop RG expansions for /3- 
functions resummed by means of generalized Pade-Borel technique result in N c3 = 3.91 |3J. 
Since the accuracy provided by this approximation was argued to be rather high (about 
1% for FP's coordinates and critical exponents), resummed e-expansion, within three-loop 
order, seems to yield plausible estimates for the quantities of interest. 



This conclusion is of prime importance. It is reasonable, therefore, to look for extra pros 
and cons it. One can find them calculating the rest of marginal values of iV from resummed 
three-loop e-expansions and making a comparison of numbers obtained with corresponding 
3D RG estimates. The inequality N c2 > 2 proven earlier |J may be used as a criterion in 
course of such a study. The e-expansions for iV cl and N c2 are found to be: 

iV c i = 2-e + ^(6C(3)-l)e 2 = 2-e+ 1.294e 2 , (5) 



iV c2 = 12-4v / 6-fl2-^) e + 
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= 2.202 - 0.569e + 0.989e 2 . (6) 
Their Pade-Borel summation gives for e = 1 

iV cl = 1.50 , iV c2 = 1.96 . (7) 

These values are close to those obtained in 3D: iV cl = 1.45, N c2 = 2.03 |J, but N c2 is 
seen to be obviously underestimated by the e-expansion since corresponding number is less 
than 2. The difference 2 — — 0.04, however, is small and may be considered as a lower 
bound for the error produced by this approximation. The most likely estimate for this error 
is believed to be close to N^ 2 D ^ — N^ 2 , i.e. being of order of 0.1. 

For iV cl the e-expansion predicts the value which is slightly greater than N^ D \ At the 
same time, these numbers differ from each other by only 3% and lie so far from the nearest 
physical value N = 2 that this difference is quite unimportant. 

We see that resummed three-loop e-expansion gives good enough numerical estimates 
for N c i and N c2 providing, however, the lower accuracy when used for evaluation of highest 
critical dimensionality N C 3. It is not surprising since the structure of series Eq. (0) turns out 
to be rather unsuitable for yielding reliable quantitative results. Indeed, to obtain precise 
numerical estimates one has to deal with series which, at least, possess coefficients decreasing 



with increasing their number. Instead, the second term in Eq. (0) exceeds, for e = 1, the first 
one. It is clear that such an expansion would not demonstrate good summability. That is 
why we believe that the true value of N C 3 is closer to 3D estimate 3.91 than to 3.39. On the 
other hand, the latter estimate differs from the former by no more than 15%. Hence, actually 
the current status of the e-expansion is not so bad provided three-loop contributions are 
properly taken into account. It seems natural that calculations of higher-order terms will 
result in further improvement of numerical estimates. 

It is worth noting that really three-loop terms added and the summation have changed 
the situation drastically. The point is that two-loop e-expansions for N c i,N c2 , and N c3 
directly extrapolated to e = 1 violate the inequalities iV cl < N c2 < N c3 which should hold 
good according to the definition of N ci ; in this approximation Nj$ < N^f < when 
e = 1. This is an alarm bell signalizing that with such short series in hand one can not 
safely penetrate into the three-dimensional world. To the contrary, the numbers given by 
the Pade and Pade-Borel resummed three-loop e-expansions at e = 1 meet abovementioned 
inequalities. 

Another point to be specially marked is as follows. Actually, we have now enough 
information resulting from higher-order RG analysis in three and 4 — e dimensions to make 
the firm conclusion about the critical behavior of the AF-like systems. Indeed, since N c2 has 
been proven to be greater than 2 |J and N c3 should exceed N c2 , magnets with planar spins 
certainly can not undergo continuous chiral transitions. Only first-order phase transitions 
into chiral states are possible, in principle, in these systems. 

We conclude with the comment concerning the ^-expansion and Monte Carlo simu- 
lations. Actually, their predictions do not contradict to those just obtained. Indeed, as 
was shown in Ref. even for simple, 0(n)-symmetric model the ^--expansion begin to 
yield reasonable numerical estimates only when n exceeds 20. Since for systems studied 
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n = 2N = 4, 6, this approach is obviously inapplicable in our case. Monte Carlo sim- 
ulations, in turn, give the values of critical exponents which are close to tricritical ones, 
especially for systems with XF-like spins [13 j. That is why it is believed 0^,0] that not 
chiral critical behavior but tricritical one or tricritical-to-critical crossover are really seen 
in these computer experiments. 

To summarize, we have found three-loop contributions to the e-expansions of critical 
order-parameter dimensionalities N c i, N c2 , and N c ^ for the model describing the chiral 
critical behavior. The Pade-Borel summation of series obtained has yielded, at e = 1, fair 
numerical estimates for iV cl and N c2 in 3D. For the lower boundary of the domain where 
continuous chiral transitions are possible the e-expansion resummed by Pade-Borel and 
simple Pade methods has given N c3 = 3.39 and N c3 = 3.81 respectively. Being close enough 
to 3D RG estimate N c ^ = 3.91, these numbers are greater than physical values N = 2 and 
N = 3. It may be considered as an evidence that in magnets with Heisenberg or planar 
spins the chiral critical behavior with specific values of critical exponents would not really 
occur and they can not belong to new, chiral class of universality. Since N c3 > N c2 and N c2 
should exceed 2, for the XY-like systems this conclusion sounds as firm. 
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